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Abstract 



Search for resonant top plus jet production in tt + jets events with the ATLAS 

detector in pp collisions at -^/s = 7 TeV 

This paper presents a search for a new heavy particle produced in association with a top or antitop 
quark. Two models in which the new heavy particle is a color singlet or a color triplet are considered, 
decaying respectively to tq or tq, leading to a resonance within the tt + jets signature. The full 
2011 ATLAS pp collision dataset from the LHC (4.7 fb _1 ) is used to search for tt events produced 
in association with jets, in which one of the W bosons from the top quarks decays leptonically and 
the other decays hadronically. The data are consistent with the Standard Model expectation, and 
a new particle with mass below 430 GeV for both W' boson and color triplet models is excluded at 
95% confidence level, assuming unit right-handed coupling. 

PACS numbers: 14.80.-j 



In the past few decades, remarkable agreement has 
been shown between measurements in particle physics 
and the predictions of the Standard Model (SM). The top 
quark sector is one important place to look for deviations 
from the SM, as the large top quark mass suggests that it 
may play a special role in electroweak symmetry break- 
ing. The recent top quark forward-backward asymmetry 
measurements from the Tevatron experiments [TJ [5] are 
in marginal agreement with SM expectations. A non- 
SM explanation could come from a possible top-flavor- 
violating process [5H5]- In these models, a new heavy 
particle R would be produced at the LHC in association 
with a top or anti-top quark. Figure [I] shows represen- 
tative production diagrams for these new particles, for 
the cases of R = W or R — <f> (see below). As shown in 
Ref. [5] , the production mechanism in pp collisions mainly 
involves quarks rather than anti-quarks at y / s=7TeV, 
even for relatively low mass particles. 

The larger number of quarks relative to anti-quarks 
produced in the initial state at the LHC leads to a res- 
onance R that decays predominantly to either the i+jet 
or t+jet final state, where baryon number conservation 
restricts the models that are available. Two models that 
can give rise to these final states are a color singlet res- 
onance (W) mostly in the tq system, and a di-quark 
color triplet model with a resonance (0) in the tq sys- 
tem. In both cases a tt+jet final state is produced, but 
a peak will be present in only one of the t+jet or i+jet 
invariant mass distributions. The new resonances are as- 
sumed not to be self-conjugate, which makes searches for 
same-sign top quarks insensitive to them [7H5]> an d to 
have only right-handed couplings. The t or i then de- 
cays to W + b or W~b, respectively. This paper considers 
the decay signature of events in which one W boson de- 
cays leptonically (to an electron or muon, plus neutrino 
final state) and the other W boson decays hadronically. 
The first direct search for such particles was performed at 
CDF [JO], which excluded color triplet resonances with 
masses below 200 GeV and W resonances with masses 
below 300 GeV, for particles with unit right-handed cou- 
pling (go) to tq. As is done in this paper, CDF used the 
formalism in Ref. [3] to define g R . CMS recently per- 
formed a search that excluded a new W' with a mass less 
than 840 GeV [TT] for particles with g R = 2 [TJ. 
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FIG. 1: Example production and decay Feynman diagrams 
for the (a) W' and (b) <f> models. 

The analysis presented here uses the full ATLAS 7 TeV 
pp collision dataset collected in 2011, corresponding to 
4.7±0.2 fb _1 of integrated luminosity [U [JJ] deliv- 
ered by the LHC. ATLAS [15] is a multi-purpose particle 
physics detector with cylindrical geometry |16j . The in- 
ner detector (ID) system consists of a high-granularity 
silicon pixel detector and a silicon micro-strip detector, 
as well as a transition radiation straw-tube tracker. The 
ID is immersed in a 2T axial magnetic field, and pro- 
vides charged particle tracking in the range \rj\ < 2.5. 
Surrounding the ID, electromagnetic calorimetry is pro- 
vided by barrel and endcap liquid-argon (LAr)/lead ac- 
cordion calorimeters and LAr/copper sampling calorime- 
ters in the forward region. Hadronic calorimetry is pro- 
vided in the barrel by a steel/scintillator tile sampling 
calorimeter, and in the endcaps and forward region by 
LAr/copper and LAr/tungsten sampling calorimeters, 
respectively. The muon spectrometer (MS) comprises 
separate trigger and high-precision tracking chambers 
measuring the deflection of muons in a magnetic field 
with a bending power of 2-8 Tm, generated by three su- 
perconducting air-core toroid systems. A three-level trig- 
ger system is used to select interesting events. The level-1 
trigger is implemented in hardware and uses a subset of 
detector information to reduce the event rate to a de- 
sign value of at most 75 kHz. This is followed by two 
software-based trigger levels, level- 2 and the event filter, 
which together reduce the event rate to ~ 300 Hz. 

Events with an electron (muon) are required to have 
passed an electron (muon) trigger with a threshold of 
transverse energy Et > 20 GeV (transverse momentum 
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Pt > 18GeV), ensuring that the trigger is fully efficient 
for the off-line selection discussed below. Electrons recon- 
structed offline are required to have a shower shape in the 
electromagnetic calorimeter consistent with expectation, 
as well as a good quality track pointing to the cluster in 
the calorimeter. Candidate electrons with Er > 25 GeV 
are required to pass the "tight" electron quality crite- 
ria |17j . to fall inside a well-instrumented region of the 
detector (\r]\ < 2.47, excluding 1.37 < \r)\ < 1.52), and to 
be isolated from other objects in the event. Muons with 
transverse momentum > 20 GeV are required to pass 
muon quality criteria |18j . to be well measured in both 
the ID and the muon spectrometer, to fall within \rj \ < 
2.5, and to be isolated from other objects in the event. 

Jets are reconstructed in the calorimeter using the anti- 
k t [12] algorithm with a radius parameter of 0.4. Jets are 
required to satisfy pt > 25 GeV and |7y| < 2.5. Events 
with jets arising from electronic noise bursts and beam 
backgrounds are rejected [20]. Jets are calibrated to the 
hadronic energy scale using p^- and ^-dependent correc- 
tions derived from simulation, as well as from test-beam 
and collision data [21]. Jets from the decay of heavy 
flavor hadrons are selected by a multivariate ^-tagging 
algorithm [55] at an operating point with 70% efficiency 
for 6-jets and a mistag rate for light quark jets of less 
than 1% in simulated ti events. Neutrinos are inferred 
from the magnitude of the missing transverse momentum 
(E™ iss ) in the event [53]. 

The signal region for this analysis is defined by re- 
quiring exactly one charged lepton and five or more jets, 
including at least one ^-tagged jet. To select events with 
a leptonically decaying W boson, events are required to 
have E% iss > 30 GeV (E^ iss > 20 GeV) in the electron 
(muon) channel. Additionally, the event must have a 
transverse mass of the leptonically-decaying W boson 
> 30 GeV in the electron channel, or scalar sum 
E^ iss + m% > 60 GeV in the muon channel [24!. Here, 
(to^) 2 = 2£ T liss £ T (l - cos<£), where E l T is the magni- 
tude of the transverse momentum of the lepton, and <fi is 
the angle between the lepton and the missing transverse 
momentum in the event. 

A variety of Monte Carlo generators are used to study 
and estimate backgrounds. The generated events are 
processed through full detector simulation [25], based 
on GEANT4 [55], and include the effect of multiple pp 
interactions per bunch crossing. To predict the event 
yield, the simulation is given an event-by-event weight 
such that the distribution of the number of pp collisions 
matches that in data. 

The ti background is modeled with MC@NLO 
v4.01 [2f interfaced to HERWIG v6.520 [53] and JIMMY 
v4.31 [29 . An additional ti sample modeled with 
MC@NLO interfaced to PYTHIA v6.425 [30] is used to 
study potential systematic uncertainties. Other ti sam- 
ples use POWHEG [3T] interfaced either to PYTHIA 
or HERWIG, as well as AcerMC v3.8 [35]. The back- 
ground from the production of single W bosons in as- 
sociation with extra jets is modeled by the ALPGEN 



v2.13 [33] generator interfaced to HERWIG. The MLM 
matching scheme [31] is used to form inclusive W bo- 
son + jets samples such that overlapping events pro- 
duced in both the hard scatter and parton showering 
are removed. In addition, the heavy flavor contribu- 
tions are reweighted using the data-driven procedures 
of Ref. [35] using the full 2011 LHC dataset. Diboson 
events are generated using HERWIG. Single-top-quark 
events are modeled by MC@NLO, interfaced with HER- 
WIG for the parton showering, in the s-channel and Wt 
channel, and by AcerMC v3.8 in the i-channel. The small 
background in which multi-jet processes are misidentificd 
as prompt leptons is modeled from a data-driven matrix 
method [36] . In determining the expected event yields, 
the ti cross section is normalized to approximate next- 
to-next-to-leading-order QCD calculations of 167+Jg pb 
for a top quark mass of 172.5 GeV [37] [35], and the total 
W+jets background is normalized to inclusive next-to- 
next-to-leading-order predictions [39] , Signal events are 
produced, for a range of W and <j) masses, with Mad- 
Graph v5.L3.16 gD] and interfaced to PYTHIA v6.425. 
Next-to-leading-order (NLO) cross sections are used for 
the predicted W boson signal normalization [6], and 
leading-order (LO) cross sections using MSTW2008 are 
used for the ^-resonance normalization [3J. 

Events are reconstructed with a kinematic fitting algo- 
rithm that utilizes knowledge of the over-constrained ti 
system to assign jets to partons. In the fit, the two top 
quark masses are each constrained at the particle level 
to 172.5 GeV by a penalty in the likelihood, computed 
from variations from this nominal value and the natural 
top quark width of 1.5 GeV. The two W boson masses 
are similarly constrained to 80.4 GeV within a width of 
2.1 GeV. This allows the z-component of the momentum 
of the neutrino from the leptonically decaying W boson 
to be computed. Both solutions from the quadratic am- 
biguity of this computation are tested when computing 
the likelihood. Charged lepton, neutrino and jet four- 
momenta are constrained in the fit by resolution trans- 
fer functions derived from simulated ti events that relate 
the measured momenta in the detector to true particle 
momenta. The full shapes of these transfer functions 
are used in the likelihood computation. All assignments 
of any four jets to partons from the ti decay are tested 
and the assignment with the largest likelihood output for 
the ti hypothesis is selected. After the assignment is se- 
lected, the originally measured jet and lepton momenta 
and E™ ss are used. The remaining jets not associated 
with the ti partons are included to form m t j and mij 
masses, where the charge of the lepton is used to infer 
which is the top candidate and which is the anti-top can- 
didate. All combinations of extra jets with the top and 
anti-top quark candidates are considered, and the pair- 
ings that give the largest mtj and mjj masses are used. 
In this way, the same extra jet can (but does not neces- 
sarily have to) be used to form mtj and mtj ■ These two 
masses are used as observables for the search. 

Several control regions are used to ensure good model- 
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ing and understanding of the backgrounds before the sig- 
nal region is examined. The preselection control region 
requires at least four jets, but does not require a b-tag. 
The dominant tt background is tested in a control region 
with exactly four jets (including at least one 6-tagged 
jet). The rejection of events with more than four jets 
reduces signal contamination. A second tt control region 
is defined by events with exactly four jets with above 
25 GeV, one of which must be 6-tagged, and exactly one 
additional jet with px between 20 GeV and 25 GeV. Signal 
contamination is further reduced by requiring that the 
AR = ^{Aif) 2 + (Acj)) 2 between the fifth jet and both 
the reconstructed top and anti-top quarks is greater than 
7r/2. Figure [2] shows distributions in the two tt control 
regions, where good agreement is observed between data 
and the prediction. The second major background, pro- 
duction of single W bosons in association with extra jets, 
is tested in a control region with five or more jets, veto- 
ing events with 6-tagged jets. The requirement of zero 
6-tagged jets reduces both signal and tt contamination. 
The distribution in Figure [3] shows good agreement be- 
tween data and the prediction within uncertainties. Ta- 
ble [I] summarizes the expected and observed yields in the 
control regions. 

Figure H] shows the expected and observed rritj and mjj 
distributions in the signal region. The data are found to 
be consistent with the SM expectation. A variety of po- 
tential systematic effects are evaluated for the predicted 
signal and the background rates and shapes. The dom- 
inant systematic effects of the jet energy scale [3T] and 
resolution [5T] lead to uncertainties of up to 10% on the 
total background rate and up to 21% on the total sig- 
nal expectation, depending on the mass of the new par- 
ticle. The other dominant systematic uncertainty from 
the difference in 6-tagging efficiency between simulation 
and data leads to uncertainties of roughly 16% on both 
the signal and background rates. Effects due to lep- 
ton trigger uncertainties and ID efficiency as well as the 
energy scale and resolution are assessed using Z — > ee 
and Z — > fifj, data, which lead to systematic uncertain- 
ties of a few percent. Other potential systematic ef- 
fects considered are the size of the small multi-jet back- 
ground (assigned 100% uncertainty); tt generator uncer- 
tainties (evaluated by comparing different results using 
the MC@NLO and POWHEG generators, 1-10%); tt 
showering and fragmentation uncertainties (evaluated by 
comparing samples using both PYTHIA and HERWIG, 
1-6%); an uncertainty on the total integrated luminosity 
(3.9%) [T51 ll4j; and the amount of QCD radiation for the 
signal and the tt background (approximately 10%, evalu- 
ated using AcerMC). Total cross section uncertainties of 
10% (55%) are used for the tt (VF+jets) backgrounds. 

Expected and observed upper limits on the signal cross 
section are computed at discrete mass points as follows. 
For each benchmark signal mass point under consider- 
ation, a signal region is defined in the rtitj-mij plane. 
When setting limits for the W' (4>) model, the m t j (mtj) 
window is significantly wider than the mtj (mtj) window 
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FIG. 2: The leading jet pt in the four-jet tt control region 
(a), and mtj in the five-jet tt control region (b). The example 
signal-only distributions are overlaid for comparison, where 
unit coupling for the new physics process is assumed. The 
total uncertainty shown on the ratio includes both statisti- 
cal and systematic effects. The "other" background category 
includes single top production, diboson production and multi- 
jet events. 
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TABLE I: Expected and observed yields in the four control regions (CR). Total refers to the total expected background, 
including ti, W+jets and the other smaller backgrounds: single top production, diboson production and multi-jet events. The 
last two lines show the expected number of events for two benchmark signal samples in each of these control regions. The errors 
include all systematic uncertainties. 





Preselection CR 


W+jets CR 


four-jet tt CR 


five-jet tt CR 


tt 


50000 ± 4700 


2000 ± 400 


19000 ± 600 


2100 ± 200 


W+jets 


46000 ± 14000 


7000 ± 2900 


3800 ± 800 


360 ± 170 


Total 


116000 ± 21000 


12000 ± 3600 


26000 ± 1300 


2900 ± 440 


Observed 


110933 


11858 


26197 


2736 


300 GeV W 


13900 ± 670 


930 ± 110 


3000 ± 400 


400 ± 80 


400 GeV (f> 


6100 ± 200 


430 ± 60 


1100 ± 100 


200 ± 20 



TABLE II: Expected and observed yields in different signal regions. The errors include all systematic uncertainties. Total refers 
to the total expected background, including ti, VK+jets and the other smaller backgrounds: single top production, diboson 
production and multi-jet events, which are not tabulated separately here. Signal window eff. refers to the efficiency for the 
signal to fall inside the optimized two-dimensional mass window. The signal region yield is calculated in the mass window at 
each benchmark signal point. Signal a refers to the total expected signal cross section, not taking into account the t (or t) plus 
jet branching fraction. 





Entire signal region 


300 GeV W 


600 GeV W 


400 GeV 4> 


800 GeV 4> 


mtj window [GeV] 




344 < m t j < 494 


566 < m tj < 904 


401 < m t j < 455 


766 < m t j < 819 


rrifj window [GeV] 




292 < m {j < 339 


549 < m { j < 650 


371 < m f j < 608 


628 < mtj < 973 


Signal window eff. 




7.5% 


9.9% 


11.9% 


5.7% 


tt 


18000 ± 3000 


740 ± 160 


270 ± 60 


660 ± 150 


60 ± 10 


W+jets 


1700 ± 560 


60 ±30 


30 ±20 


80 ±40 


8±5 


Total 


22000 ± 3700 


820 ± 190 


320 ± 80 


780 ± 180 


70 ±20 


Observed 


22731 


970 


343 


923 


77 


Signal region yield 




560 ± 120 


98 ± 24 


410 ± 100 


20 ±6 


Signal a 




19.0 pb 


1.55 pb 


7.9 pb 


0.67 pb 



to account for the fact that the resonance is predomi- 
nantly in the mjj (mtj) system. The windows are op- 
timized to maximize sensitivity, accounting for the full 
effect of systematic uncertainties. Typical mass win- 
dows are shown in Table [n] For each mass window, 95% 
confidence level (C.L.) upper limits on the signal cross 
section (times the branching ratio to t or t) are com- 
puted using a single bin frequentist CL S method [4"2] . 
No shape information is used within the mass windows. 
Table [TT] shows the expected and observed event yields 
in several of the signal region windows. Expected and 
observed 95% C.L. lower limits on the signal mass are 
derived, assuming a coupling of — 1 and gu — 2, and 
are shown in Figure [5] Assuming that the cross section 
scales as g\, the exclusion in the mass-coupling plane is 
shown in Figure [6] As shown, most of the parameter 
space in this model, which was favored by the Tevatron 
forward-backward asymmetry and cross section measure- 
ments [43], has been excluded. 

In conclusion, this paper presents a search for a new 
heavy particle R in the tj or tj system of ti plus extra 
jet events with the ATLAS detector. Such new parti- 
cles have been proposed as a potential explanation of the 
difference from the SM values of the forward-backward 
asymmetries measured in top quark pair production at 
the Tevatron. The full 2011 ATLAS pp dataset (4.7 fir 1 ) 
is used in the search. Assuming unit coupling, the ex- 
pected 95% C.L. lower limit on the mass of the new par- 



ticle is 500 (700) GeV in the W (</>) model. No signif- 
icant excess of data above SM expectation is observed, 
and 95% C.L. lower limits of 430 GeV for both the W 
and 4> models are set. At = 2, the limits are 1.10 
(1.45) TeV for the W (</>) model, with expected limits 
of 0.93 (1.30) TeV. These are the most stringent limits 
to date on such models. Most of the regions of parame- 
ter space for these models that are more consistent with 
the Tevatron forward-backward asymmetry and ti cross 
section measurements than the SM are excluded at 95% 
C.L. by these results. 
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FIG. 3: Expected and observed distribution of m t j in the 
W+jets control region. The example signal-only distributions 
are overlaid for comparison, where unit coupling for the new 
physics process is assumed. The total uncertainty shown on 
the ratio includes both statistical and systematic effects. The 
"other" background category includes single top production, 
diboson production and multi-jet events. 
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FIG. 4: Expected and observed distributions of (a) m t j and 
(b) mtj in the signal region. The example signal distribu- 
tions assume unit coupling for the new physics process. The 
total uncertainty shown on the ratio includes both statisti- 
cal and systematic effects. The "other" background category 
includes single top production, diboson production and multi- 
jet events. 
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FIG. 5: Expected and observed 95% C.L. upper limits on 
the (a) W' and (b) <j> model cross sections. The CDF result 
is documented in Ref. [10] . The W' cross sections are NLO 
calculations, and the <f> cross sections are LO calculations. 



FIG. 6: Expected and observed 95% C.L. upper limits on 
the (a) W 1 and (b) 4> model cross sections assuming a cross 
section which scales with g\. The hatched area shows the 
region of parameter space excluded by this search at 95% 
C.L. The CDF result is documented in Ref. [10]. The W 
cross sections are NLO calculations, and the (j> cross sections 
are LO calculations. The region favored by the Tevatron Afb 
and a t i measurements is shown as the dark band [43 . 
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